Gravitational Optics: 
Self-phase modulation and harmonic cascades 

J.T.Mendongag 

GoLP, Instituto Superior Tecnico, 
1049- 001 Lishoa, Portugal 

^ \ Vitor Cardoso[| 

§ ■ CENTRA, Departamento de Fisica, Instituto Superior Tecnico, 

> 
O 

^ : Abstract 



O 
u 



Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 
(Dated: February 3, 2008) 



Nonlinear wave interaction of low amplitude gravitational waves in flat space-time is considered. 
Analogy with optics is established. It is shown that the flat metric space-time is equivalent to a 
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Q\ • centro-symmetric optical medium, with no second order susceptibility. The lowest order nonlinear 

o . 

CN , effects are those due to the third order nonlinearity and include self-phase modulation and high 



harmonic generation. These processes lead to an efficient energy dilution of the gravitational wave 
energy over an increasingly large spectral range. 
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I. INTRODUCTION 



Gravitational radiation is a direct result of the theory of gravitation |l|, |] and, even if not 
yet directly observed, its existence has been indirectly inferred from the change of orbital 
parameters in binary star systems 0. Potential gravitational wave sources 0] can emit 
intense and short bursts, such as collapsing massive binary stars or supernova explosions, and 
continuous waves of much lower amplitudes, such as stable binary stars or non-axisymmetric 
spinning stars. 

It is well known that the Einstein's equation governing the gravitational field contains 
strong nonlinearities. The nonlinear effects will, of course, be more relevant to the short 
bursts, containing only a few cycles of radiation, and associated with the more catastrophic 
events. It is also known that nonlinear gravitational waves of the soliton type can be excited 
@) H) 0) Hi- We have learned, however, from the more common situations of nonlinear 
optics [^, that the excitation of solitons depends very critically on the balance between 
the dispersive and the nonlinear effects. Even if the gravitational solitons do not seem 



to display such a balance and behave much more like superpositions of linear waves [10 



they nevertheless correspond to particular solutions. And, with the exception of the soliton 



representation of cosmological solutions |TT|, |T2[ |T3[, which are pertinent to the large scale 
structure of the universe, it is very unlikely that solitons will spontaneously be excited by 
local gravitational wave sources. 

It is therefore quite plausible in physical terms to consider the problem of gravitational 
wave-packets of arbitrary shape. This will be considered here. In order to understand their 
typical behavior, and to describe the main physical pictures of their spectral evolution, we 
concentrate here on the simplest possible situation: that of propagation in a flat space-time. 
We are not interested in the wave generation phenomena, for which there exists a standard 



powerful formalism |]I4[ valid even in strongly curved regions, nor in the propagation of 
gravitational waves in strongly curved regions. Here we want to focus on regions far from 
the sources, where the background is approximately fiat. This allows us to simplify consid- 
erably the problem but still retaining some important physical features. We will try then 
to establish clear analogies with the well known concepts of nonlinear optics. In particular, 
the nonlinear gravitational susceptibilities will be established and the nonlinear wave mix- 
ing processes will be considered. Explicit analytical solutions for the processes of self-phase 
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modulation and of high-order harmonic cascades will be derived. 



II. BASIC EQUATIONS 

Einstein's equation for the gravitational field in the absence of matter is given by: 

p -pi -pi I p^ -pm-pl rv /-I \ 

^ik — ^ ik,l ^il,k^^ik^lm ^ il km — K^J 

where the comma stands for ordinary derivative, and the Christofell symbols are determined 
by the derivatives of the metric tensor elements, according to: 

We are considering a flat space-time, perturbed by a small amplitude gravitational wave. 
This means that we can write: 

9ij ~ Vij ~l~ ^iji (3) 

where \hij\ <IC 1, and: 

^700 = 1 , Vii = -1 {i = 1, 2, 3) , rjij = {i^ j). (4) 
In this case, we can write: 

Rik = Rik + Rik^: (5) 

where the linear term is: 

j^L _ ^ Im f _|_ d'^hmk _ d'^hik _ d'^hml \ /gs 

2^ ydx'-dx^ dx^dx'^ dxWx"^ dx'^dx^ J 
The nonlinear term Rf^ contains second, third and higher order nonlinearities. This 
means that, in principle, we can have three, four and higher wave mixing processes. Let us 
then write: 

R^ = R^ + R^ + -, (7) 

where we have, for the lowest order term: 
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The third order nonhnearities are contained in: 
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We could go on to higher orders but as we assume \hij\ ^ 1, these terms become less 
important. Furthermore, stopping at third order suffices, as we will show, to uncover some 
truly non-linear and interesting aspects of gravitation. 

It is well known that the linear term (P) can simply be written as: 



Kk ^ -i^d'd.K,, 



(10) 



where the d'Alembert operator is: 



d^d,, 
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dx^dxj 



n 



d^ 



dxWx^ 

This means that we can write the nonlinear wave equation as: 



(11) 



d^d„h 



'p'Hk 



p(2) , p(3) , d(4) 
^ik +^ik +^ik 



(12) 



III. NONLINEAR WAVE COUPLING 



Equation ([1^) predicts the possibility of several types of nonlinear gravitational wave 
coupling. This can be studied by assuming that several waves of the type: 



hij{n) = eij{n)a{n) exp[iqk{n)x^], (13) 

coexist in the same region of flat space-time. Here a{n) are slowly varying amplitudes, 
eij(n) are unit polarization tensors such that e*jeij = 1, and qk{n) are the four-wavevectors 
associated with the interacting waves n = 1, 2, 3... 

The total wave field, for real waves, will be given by: 
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hij = ^ hij{n) + c.c. (14) 

n 

In order to determine the evolution of the amphtudes a{n), we can construct envelope 
equation, by noting that: 

d^d^hij{n) ~ (^qtq^ + ^^k-^^ hij{n). (15) 
Assuming that the waves still obey the linear dispersion relation, valid in the limit Rfj^ 0: 

qkin)q\n) = 0, (16) 

we obtain, for three wave mixing: 

2gfe(l)^a(l) = vil)a{2)a{3) expitAq^x"), (17) 



where: 



Agfc = gfc(2)+gfe(3)-gfc(l). (18) 

The nonlinear coupling coefficient v{l) is determined by Rif, equation(H). It can be 
verified that v{l) = 0, which means that the three wave coupling is forbidden. We can 

(2) 

use the hard way, by calculating explicitly R^j , or the easy way, by noting that fiat space- 
time is a centro-symmetric medium and therefore, like in optics, it has no second order 
susceptibility. We then proceed to the next order of nonlinearity, and write the envelope 
equation for the four wave mixing process: 

zgfe(l)^a(l) = w{l)a*{2)a{3)a{A) expizAq^x''), (19) 

where now: 



Agfc = qk{3) + gfc(4) - gfe(l) - ^^(2), (20) 

and the coupling coefficient ti'(l) is determined by Rif ■ Similar equations given the other 
three wave amplitudes. In order to derive explicit results, we assume wave propagation along 
a given direction x^ = x. We have: iqkx'' = iq^x^ + iqiX^ . We will use: qi = q, qo = — fi/c, 
and x^ = ct, which then leads to: 
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exjp(iqkX^) = exp{iqx — i^lt). (21) 
The hnear dispersion relation (|TB|) is: 



Qkx'^ = ^ - g2 = 0. (22) 

This shows that flat space-time is non-dispersive. This means that, for Q{2) + f2(l) = 
f2(l) + fi(2), we can guarantee a perfect phase matching: 

Agfc = 0. (23) 

We will also have: 

dxk ^^'^^ dx^ ^ \dx ~^ c dt ) ^ '^^ 

We can now introduce a variable z = x — ct, such that: 

d 

iq{l)—a{l) = w(l)a*(2)a(3)a(4). (25) 
Let us use equation (|) to calculate. The result is: 



w{l) = ie*fc*,^«-e'P*(2) {[gK3)e™(3) + q^{3)e^l{3) - qU^)e,i{3)] ■ 
[g„(4)epfc(4) + gfc(4)ep„(4) - gp(4)efe„(4)] - gi(4)enp(4) 
fefe(3)epi(3) + gi(3)epfc(3) - gfc(3)eifc(3)]} + permMtatzon(2, 3, 4). (26) 

This is a very complicated expression, where the last term represents two terms formally 
identical to the first one, and obtained by permutation of the indices 2, 3 and 4. It is 
important to notice that, in order of magnitude, we have w{l) ~ Vt^ 



IV. SELF-PHASE MODULATION 



An important particular case is obtained when the modes n = 1,2,3,4 coincide. The 
mode coupling equation (|^) is then reduced to: 

d 

iq-^T-a = wlapa. (27) 
oz 
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It is instructive to compare it with the envelope equation for an optical pulse (for instance 
a laser pulse), with central frequency c^o, propagating in a nonlinear medium with group 
velocity Vq. Its electric field amplitude is described by the equation 0: 

d 1 d\ 

— + ——] Eo = iuJoa\Eo\'^Eo. (28) 



^ dx vq dt 

where Eq is the electric field amplitude and a is the nonlinear coupling coefficient 



Notice that equation can be written in a similar form: 



H — — a = —i—\a\ a. (30) 



dx cdt J q 

Comparing with equation (|2^) we conclude that the equivalent nonlinear susceptibihty 
of the fiat space-time, which would be the gravitational wave version of the electromagnetic 
susceptibility x^^"* is equal one, whereas in optical materials the typical value is x*-^"* ~ 10~^^. 
This means that he gravitational field is extremely nonlinear and that we expect to observe 
the same nonlinear effects as in optics at a much lower wave amplitudes. In particular, we 
can follow the usual optical approach and obtain a wavepacket solution of the form: 

a{x,t) = a{z) exp [i(f){z,t)] , (31) 

where the nonlinear phase is: 

TV 

(j){z,t) = (f)o--\a{z)\H. (32) 

This means that a gravitational wavepacket, even of very small amplitude, will suffer a self- 
phase modulation such that its central frequency Q will not remain fixed but will change 
along propagation, according to: 

n = n{o) - ^ = n{o) - -t|-|a(^)p. (33) 

at qc oz 

For instance, if we have a gaussian gravitational wave pulse of the form: 

|a(^)r = «oexp f-^V (34) 
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we will obtain a gravitational wave frequency shift AQ = Q — Q{0): 



An = -'^\a{z)\H. (35) 
qc 

It contrast with what usually occurs in nonlinear optics, this frequency shift will be 
positive at the pulse front (2; > 0), and negative at the pulse rear {z < 0). We also see that, 
for very large propagation times t this shift can be quite significant, AQ ^ ^(0), even for 
weak gravitational perturbations, a <^ 1, due to the already noticed strong nonlinearities. 
To be quantitative, let's see what conditions have to be fulfilled in order to have the same 



An An 

Hit; bfciiiit; litjqutjiii;^ smii 

to satisfy this condition we have to have 



effect as in nonlinear optics, i.e., to have the same frequency shift — = "grav" order 



|a(.)r~|Eo|\S;^, (36) 

"'grav 

where dopt is a length travelled by light in usual experiments, and dgrav is the length travelled 
by a gravitational wave until it reaches the gravitational wave detector. As a typical example, 
we have dopt ~ lO^^m, |-EoPXopt ~ 10~^. For these values, we get 

6 X 10-^=^ 

[dgrav{KpC)\ 

where the distance is now measured in Kpc. For a gravitational wave produced by any 
reasonable astrophysical event, our best expectations |jT^ give a ~ 10~^^, as the wave arrives 



on Earth. But, along its path from the radiation source to the Earth, its amplitude will 
remain significantly higher over large distances. It seems therefore very likely that such an 
effect will indeed take place. 



V. HARMONIC CASCADE 



Another possibility of nonlinear gravitational wave coupling is the generation of a large 
spectrum of harmonics of the initial frequency Q{0). Here, in principle, the sum over n will 
extend to infinity: 

n=— 00 

We will also have: 
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Qkin) = nqk. (39) 

Four wave mixing between these different harmonics will occur, where a perfect phase- 
matching condition can be verified: 



Agfe = [qk{2) + qk{n - 1)] - [^^(1) + qk{n)] = 0. (40) 

For the present case, an envelope equation for the amplitudes of the different harmonics 
can be derived, in the same way as equation ([l9|). Assuming that the fundamental and the 
first harmonics are dominant, \hij{n = 1,2) ^ hij{n ^ 1,2), we can retain only two of the 
mixing terms and get, for propagation along one given direction x: 

\ dx^~cdt) ^^^^ ^ ^^^"^ [a*{l)a{2)a{n - 1) + a{l)a*{2)a{n + 1)] . (41) 



where w{n) is determined by an expression similar to equation (|26|) . Let us now define: 



a*(l)a(2) = Ie'\ (42) 

where / and 6 are real. In order to get order of magnitude analytical solutions, let us assume 
that w{n) oc q{n), and let us also neglect the amplitude variations of the dominant modes, 
which means that the intensity parameter I is approximately constant. We are then led to: 

d 

— a(n) = iw \e'^a(n - 1) + e-'^a(n + 1)1 , (43) 
oz '- -' 

with: 

win) 

w = ^J. (44 

q{n) 

It is now convenient to introduce new amplitude variables b{n), such that: 

a{n) = (- 1)^^/2 gm<5^^^)^ (45) 

and we get: 

■^ = w [b{n - 1) - b{n + 1)] . (46) 
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Using T = wz, we can reduce this equation to the well known recurrence relation for Bessel 
functions: 

^b{n) = b{n-l)-b{n + l). (47) 

An adequate solution will then be given in terms of Bessel functions of the first kind 
Jre(T), such as: 

b{n,T)=AJ4T)+A'Jn-l{T), (48) 

where the two constants A an A' are determined by the initial conditions: 

A^ a{n^l,T = 0) , A' ^a{n = 2,T = 0), (49) 

and all the other harmonic wave amplitudes as assumed to be initially equal to zero: a{n ^ 
1, 2;t = 0). Notice that this solution satisfies energy conservation, in the sense that: 

oo 

\b{n,r)\'^ = const. (50) 

n=—oo 

It is then an adequate, if if particular, solution to the problem of harmonic generation of 
low amplitude gravitational waves propagating in a flat space-time. 



VI. CONCLUSIONS 



Nonlinear wave coupling of gravitational waves was considered in this work. The problem 
of low amplitude waves in flat space-time was examined, in order to extract the main physical 
consequences with the simplest possible formal complexity and to establish clear connections 
with the well known effects occurring in nonlinear optics. 

In particular, the gravitational equivalent to the optical second and third order suscepti- 
bility was derived. It was shown that the second order gravitational susceptibility is equal 
to zero and, consequently, the process of three nonlinear wave coupling is forbidden. This 
is not surprising because the flat space-time background considered here can be seen as 
a centro-symmetric material medium. This property will eventually disappear in curved 
space-time. 
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In contrast, the gravitational third order susceptibihty was shown to be equal to one, 
thus meaning that the fiat space-time is a strongly nonlinear medium. As in optics, four 
wave mixing processes can then be considered. But, unlike the case of optical phenomena, 
these wave mixing processes can become relevant even for very low wave amplitudes. The 
particularly important cases of generation of harmonic cascades and of self phase modulation 
were considered, and explicit analytical solutions were derived. They show that a significant 
spectral broadening and spectral energy dilution can then take place, specially in the vicinity 
of radiation sources where the wave amplitudes are non- negligible. Such energy dilution is 



not considered in the usual estimates for gravitational wave detection 0, |15[ and could 
eventually lower the prospects for direct observation of these waves near the Earth. 

Finally, it should be noticed that self phase modulation occurring in optics can be de- 



scribed as a particular example of photon acceleration |16, 17]. In analogy, the present result 
showing the existence of self phase modulation of gravitational wave pulses could also be 
interpreted as graviton acceleration. This aspect will be explored in a future work. 
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